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ABSTRACT 

The  XFV-12A  incorporates  thrust  augmentation  to  generate 
propulsive  lift.   Under  the  influence  of  forward  flight  the 
lift  jets  are  subject  to  bending,  which  produces  shed  vor- 
tices that  induce  velocities  at  the  aircraft.   Additionally, 
operations  in  ground  effect  during  hover  and  for  small  for- 
ward velocities  can  result  in  reingestion  of  engine  exhaust 
gases.   The  objectives  of  this  thesis  were  twofold.   The 
first  was  to  investigate  analytically  aspects  of  the  vertex 
shedding  phenomenon  through  the  application  of  potential  flow 
around  the  lifting  surfaces  of  the  aircraft.   Secondly,  ex- 
perimental analysis  was  performed  to  determine  the  unique 
reingestion  properties  of  the  XFV-12A  operating  statically 
under  various  wing  and  canard  configurations. 

The  pressure  coefficient  for  jet  curvature  obtained  under 
two-dimensional  potential  flow  was  8.66.   Reingestion  of  ex- 
haust gases  was  primarily  through  the  main  inlet  and  not  the 
auxiliary  inlet;  the  source  of  most  of  the  exhaust  gas  was 
the  canard. 
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I .   INTRODUCTION 

The  XFV-12A  is  a  V/STOL  aircraft  currently  under  develop- 
ment by  the  Columbus  Aircraft  Division  of  North  American 
Rockwell.   Its  design  incorporates  high  performance  super- 
sonic operation  in  a  fighter  or  attack  mission  as  well  as 
vertical  and  short  take-off  operations  from  proposed  ships 
of  the  1980's.   A  unique  property  of  the  XFV-12A  is  incor- 
poration of  thrust  augmentation  when  operating  under  propul- 
sive lift.   Combining  the  effects  of  hypermixing  nozzles  and 
Coanda  blowing,  this  aircraft  is  designed  to  take  off  ver- 
tically with  a  thrust-to-weight  ratio  of  less  than  one. 
When  it  becomes  operational,  the  XFV-12A  will  offer  the  first 
employment  of  thrust  augmentation  in  military  or  civilian 
aircraft . 

Problem  areas  under  investigation  include  the  shedding 
of  vortices  from  the  exhaust  jets  and  an  experimental  analysis 
of  the  XFV-12A's  recirculation  properties.   The  vortex  shed- 
ding phenomenon  occurs  when  a  downflow  of  the  lift  jets  is 
modified  by  forward  velocity,  ground  effect,  or  both.   The 
resultant  vortices  are  shed  as  a  result  of  a  bending  of  the 
lift  jets  as  shown  hypothetically  in  Figure  1.   This  bending 
is  derived  from  a  high  pressure  area  underneath  and  a  low 
pressure  area  above  the  jet.   These  vortices  are  generated 
in  much  the  same  manner  as  wing  tip  vortices  produced  in 
finite  airfoil  theory.   The  nature  and  influence  of  these 
vortices  are  the  unknowns  under  investigation. 
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A  potential  flow  solution  to  a  propulsive  lift  system 
was  performed  by  Albers  and  Potter  in  1971  [Ref.  1] .   This 
analysis  assumed  an  exhaust  profile  in  accordance  with 
Spence's  theory  [Ref.  2]  for  a  jet  augmented  flap.   The 
McDonnell  Douglas  Corporation  [Ref.  3]  has  conducted  work 
on  a  jet-wing  lifting  surface  theory  using  elementary  vor- 
tex distributions.   A  solution  for  the  exhaust  profile  it- 
self, however,  is  not  known. 

The  vortex  shedding  is  a  consequence  of  the  pressure 
difference  on  the  two  sides  of  a  bent  jet.   To  gain  insight 
into  the  vortex  problem,  a  potential  flow  model  was  formu- 
lated to  determine  jet  bending.   The  potential  flow  was 
solved  by  two  methods.   A  computer  program  was  used  to  solve 
Laplace's  equation  for  two-dimensional  potential  flow.   For 
this  case,  only  the  canard's  influence  was  considered.   A 
continuous  solution  was  obtained  using  an  electrostatic 
analogy.   Conducting  paper  provided  the  flow  domain.   Laplace's 
equation  was  again  applied,  but  under  the  influence  of  both 
the  wing  and  canard  in  a  two-dimensional  analysis. 

Recirculation  experiments  provide  valuable  data  which 
adequately  depict  the  reingestion  properties  of  the  XFV-12A. 
Since  this  characteristic  is  not  readily  calculated,  results 
are  based  purely  on  experimental  work  performed  and  the  cur- 
rent literature  available  on  the  subject.   In  subsequent 
discussions,  the  terms  reingestion  and  recirculation  will 
be  used  to  describe  the  same  phenomenon. 

Considerable  work  has  been  performed  in  the  area  of 
recirculation  properties  of  lift  engines.   Investigations 
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of  recirculation  for  lift  engines  of  various  configurations 
applied  to  various  airframe  shapes  have  been  performed 
[Refs.  4,5].   An  analysis  to  determine  the  applicability  of 
small  scale  testing  has  been  done  [Ref.  6].   Reduction  of 
reingestion  was  examined  by  applying  several  concepts  to 
prevent  recirculated  gases  from  entering  lift  engine  inlets 
[Ref.  7] .   Conclusions  have  not  always  been  consistent  when 
analyzing  the  most  significant  contribution  to  recirculation 
or  the  optimal  method  of  prevention.   However,  it  was  commonly 
held  that  reingestion  is  a  unique  property  for  each  airframe 
and  accompanying  engine  configuration.   Thus  the  particular 
properties  of  the  XFV-12A  operating  in  the  vertical  or  short 
take-off  mode  offer  a  unique  problem. 

The  reingestion  problem  was  investigated  experimentally. 
The  literature  was  helpful  in  formulating  and  designing  ex- 
perimental apparatus  and  procedures.   For  carrying  out  the 
experiments,  a  one -tenth  scale  semi -span  model  of  the  XFV-12A 
was  used.   In  order  to  simulate  a  zero  forward  velocity, 
testing  was  performed  statically.   Recirculation  properties 
of  interest  included  results  for  both  the  STOL  and  VTOL 
configurations  at  varying  heights  above  the  ground.   In  all 
cases,  the  reingestion  at  both  the  main  and  auxiliary  inlets 
was  determined. 

The  computer  program  and  its  background  are  discussed 
in  Section  II  of  this  thesis.   Section  III  contains  the  work 
done  with  the  electrostatic  analogy.   The  next  section  describes 
the  arrangement  of  the  experimental  apparatus  including  supple- 
mental hardware  necessary  to  supply  appropriate  simulation  of 
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the  lift  engines.   Section  V   covers  experimental  work  per- 
formed preliminary  to  the  actual  reingestion  tests.   The 
reingestion  experiments  are  described  in  Section  VI.   Ex- 
perimental results,  conclusions,  and  recommendations  follow. 
Graphical  test  results  and  test  data  are  recorded  subsequent 
to  the  text. 
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I I .   COMPUTER  SOLUTION  TO  THE  POTE NTIAL  F LOW  EQUATIONS 
FOR  TWO-DIMENSIONAL  THRUST  AUGMENTOR" CANARD 


A.   ASSUMPTIONS 

In  formulating  the  problem,  the  mass  flow  rate  of  the 
air  leaving  the  diffuser  exit  plane  was  assumed  to  be  equal 
to  the  secondary  air  entering  the  diffuser.   This  assumption 
was  essential  for  satisfying  continuity  in  the  control  volume 
which  encircled  the  canard.   (Using  an  experimentally  ob- 
tained augmentation  ratio  of  6.84  [Ref.  8],  the  actual  ratio 
of  exhaust  at  the  canard  exit  to  the  secondary  flow  is  about 
eight  to  seven.) 

The  working  fluid  was  assumed  to  be  incompressible,  ir- 
rotational,  and  two-dimensional.   With  the  irrotationality 
assumption,  Euler's  equations  of  motion  for  fluids  become 
homogeneous.   When  a  single  dependent  variable  representing 
the  potential  function  is  introduced,  a  single  partial  dif- 
ferential equation  results. 

8x2    3y2  v  J 

The  conjugate  function  to  the  potential  function,  or  the 
stream  function,  is  derived  from  the  continuity  equation  for 
incompressible  flow. 

|~f  +  pL   =  0.  (2) 

8x2    8y2  v  ' 

The  equation  for  the  stream  function  was  solved  in  the  com- 
puter program. 
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B.   NATURE  OF  PROBLEM 

The  equation  for  the  stream  function  is  elliptical  in 
nature  yielding  a  boundary  value  problem.   The  boundary 
conditions  were  of  the  Dirichlet  type.   Values  of  the  stream 
function  were  specified  on  all  edges  of  the  domain.   Values 
of  the  stream  function  were  also  assigned  to  the  canard 
flaps  [Ref .  9] .   The  boundary  conditions  were  held  constant 
on  the  upper  and  lower  bounds.   The  rate  of  change  of  stream 
function  values  on  the  upstream  and  downstream  boundaries 
was  held  constant.   These  boundary  conditions  duplicated 
freestream  conditions. 

For  simplicity,  only  the  canard  was  considered  in  this 
analysis.   It  was  positioned  nine  chord  lengths  from  the 
front  boundary  of  the  domain  to  its  leading  edge.   The  rear 
of  the  domain  was  positioned  about  40  chord  lengths  aft  of 
the  canard.   At  this  position  it  was  anticipated  that  the 
flow  field  disturbances  generated  by  the  canard  would  have 
returned  to  freestream  conditions  at  the  boundaries.   The 
discretized  domain  is  illustrated  in  Figure  2. 
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C.   DISCRETIZATION 

Since  the  canard  is  small  compared  to  the  size  of  the 
domain,  a  variable  mesh  was  necessary  to  represent  optimally 
the  flow  while  economizing  on  computer  storage  requirements. 
It  was  anticipated  that  the  direction  of  the  streamlines 
would  be  changing  most  rapidly  in  the  vicinity  of  the  canard, 
The  highest  concentration  of  solution  points  was  then  near 
the  canard.   The  coarsest  distribution  of  solution  points 
was  located  near  the  corner  most  distant  from  the  canard. 
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The  final  distribution  of  solution  points  within  the  domain 
resulted  in  201  distinct  points. 

Using  a  Taylor  series  expansion,  the  differential  equa- 
tion can  be  reduced  to  a  difference  equation  as  illustrated 
in  Appendix  A.   Each  solution  point  is  expressed  in  terms 
of  the  four  points  adjacent  to  it.   The  resultant  error  on 
transition  to  difference  equations  is  of  the  fourth  order. 
Deriving  a  difference  equation  for  each  solution  point  re- 
sults in  201  simultaneous,  linear,  algebraic  equations  in 
201  unknowns.   In  matrix  form,  the  problem  is  easily  stated 
as 

A-x  =  B  (3) 

where  A  is  the  matrix  of  coefficients  for  the  unknown  values 
of  the  stream  function,  x  is  the  vector  of  unknown  solution 
points,  and  B  is  a  vector  containing  the  boundary  conditions 
assigned  to  both  the  perimeter  of  the  domain  and  the  canard 
flaps . 

Since  most  of  the  elements  in  A  are  zero,  it  became  more 
efficient  to  utilize  a  subroutine  requiring  only  the  storage 
of  the  nonzero  elements  and  as  few  zero  elements  as  possible 
Such  a  library  subroutine  in  the  IBM  360  system  at  the  Naval 
Postgraduate  School  is  "GELB."   This  subroutine  employs 
Gaussian  elimination  techniques  in  solving  a  banded  matrix. 
This  reduced  the  storage  requirements  for  this  particular 
problem  by  a  factor  of  about  ten.   Not  only  did  "GELB"  re- 
duce the  time  required  to  obtain  a  solution,  but  it  provided 
for  a  finer  mesh  size  within  the  domain  for  a  given  quantity 

of  core  space . 
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III.  ELECTROSTATIC  ANALOGY  SOLUTION  TO  THE  POTENTIAL  FLOW 
EQUATION  FOR  A  TWO-DIMENSIONAL  THRUST  AUGMENTOR  WING 
AND  CANARD  SYSTEM 


A.   ELECTROSTATIC  ANALOGY 

The  differential  form  of  Gauss'  flux  theorem  as  it 
governs  electric  fields  states  that  the  divergence  of  the 
field  is  proportional  to  the  charge  density  [Ref .  10] .   In 
vector  notation: 

V-E  =  p/eo.  (4) 

Since  the  electric  field,  E,  at  any  point  is  equal  to  the 
negative  of  the  gradient  of  the  electric  potential,  V, 

E  =  -  VV.  (5) 

Gauss'  flux  theorem  can  be  rewritten  as  Poisson's  equation. 

V2V  =  "P/eo.  (6) 

If  Poisson's  equation  is  assumed  to  apply  in  a  region 
of  zero  charge  density,  the  expanded  two-dimensional  electric 
potential  equation  can  be  written: 

This  equation  is  Laplace's  equation  as  it  applies  to  elec- 
tric fields  [Ref.  10] . 

If  the  voltage  is  allowed  to  be  proportional  to  the 
stream  function,  a  plot  of  equipotential  lines  on  some  con- 
ducting medium  will  yield  a  continuous  solution  for  the 
stream  function  [Ref.  11].   Additionally,  the  zero  charge 
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density  assumption  is  appropriate  since  it  is  analogous  to 
the  irrotationality  requirement  of  potential  flow. 

Similarly,  the  current  flowing  in  the  domain  is 
analogous  to  the  fluid  potential  function.   Since  circula- 
tion is  related  to  the  potential  function  by 

r  =  /  d<f>  (8) 

c 

the  value  of  the  current  can  be  used  to  measure  circulation 
[Ref.  11]. 

B.   EXPERIMENTAL  APPARATUS 

The  hardware  necessary  to  provide  such  a  solution  in- 
cluded conducting  paper  of  a  known  resistance,  an  electro- 
meter for  use  in  measuring  voltages  and  currents,  a  ten-volt 
d-c  power  source,  and  seven  potentiometers  (Figure  3).   The 
potentiometers  took  a  predesignated  fraction  of  the  ten-volt 
source  and  applied  it  to  a  particular  component  of  the  wing- 
canard  system.   A  special  pen  which  contained  a  conducting 
paint  was  used  to  draw  the  wing  and  canard  onto  the  conduct- 
ing paper. 

The  resistance  of  the  paper  was  determined  using  a  known 
voltage  source  and  the  electrometer  which  measured  current. 
The  electrometer  was  placed  in  series  with  the  conducting 
paper  and  the  voltage  source.   The  resistance  was  measured 
to  be  approximately  10,000  ohms  per  centimeter.   The  poten- 
tiometers were  selected  so  that  their  resistances  would  be 
much  less  than  that  of  the  paper,  but  large  enough  to  keep 
a  low  amount  of  current  flowing  from  the  ten- volt  source. 
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If  the  maximum  allowable  current  in  the  battery  was  con- 
strained to  be  ten  milliamps,  then  the  potentiometers  had. 
to  be  on  the  order  of  1,000  ohms.   This  satisfied  both  the 
current  and  resistance  constraints. 

The  voltage  induced  on  the  wing  and  canard  flaps  were 
assigned  so  as  to  provide  for  equal  circulation  around  the 
forward  and  aft  flaps  of  the  wing  and  canard.   In  addition, 
voltages  were  constrained  to  provide  a  flow  velocity  in  the 
canard  and  wing  augmentors  which  was  twice  that  of  freestream 
conditions . 

The  equal  circulation  requirement  was  implemented  by 
providing  equal  but  opposite  current  flowing  in  the  two 
Coanda  flaps.   The  flow  velocity  constraint  was  imposed  by 
comparing  the  distance  between  the  Coanda  flaps  on  the  paper 
with  the  distance  between  the  streamlines  upstream  which 
stagnated  on  those  Coanda  flaps. 

The  center  ejector  flaps  were  constrained  to  have  zero 
current  and  a  voltage  value  midway  between  those  assigned 
to  the  Coanda  flaps.   These  requirements  resulted  in  a 
negligible  lift  contribution  made  by  the  center  ejectors  as 
well  as  an  equal  distribution  of  wing  and  canard  augmentor 
inlet  flow  on  each  side  of  the  center  ejector. 

In  conducting  the  experiment,  however,  the  equal-but- 
opposite  circulation  requirement  could  be  imposed  on  the 
canard  only.   The  results  shown  in  Figure  4  indicate  the 
streamlines  having  the  same  voltage  value  as  those  induced 
on  the  wing  and  canard  flaps.   The  ground  is  the  only  por- 
tion of  the  boundary  depicted. 
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I V .   EXPERIMENTAL  DESIGN  AND  MODEL  MODIFICATION 

A.   PLUMBING  AND  FIXTURES 

Compressed  air  was  supplied  by  two  compressors  rated  at 
17.2  and  49.5  cubic  feet  per  minute,  respectively.   One  was 
equipped  with  a  100  gallon  storage  tank,  while  the  other  had 
an  80  gallon  tank.   These  compressors  were  installed  in  paral 
lei  with  three  additional  storage  tanks  with  a  capacity  of 
400  gallons  each.   The  total  storage  capacity  of  the  system 
was  1,380  gallons  or  181.3  cubic  feet. 

The  compressors  were  limited  by  design  to  a  pressure  of 
137  lb/in2.   The  actual  operating  pressure  of  the  system  was 
constrained  to  be  no  more  than  125  lb/in2.   For  the  planned 
operating  velocity  at  the  ejectors  of  136  ft/sec,  substantial 
ly  more  than  15  minutes  of  experimental  run  time  at  the  de- 
sired operating  conditions  was  available. 

A  single  regulator  and  stop  valve  connected  the  storage 
tanks  to  the  plenum,  which  supplied  the  ejector  hypermixing 
nozzles  and  Coanda  blowing.  This  provided  for  single-valve 
control  of  the  experiments.  Once  the  desired  pressures  had 
been  obtained  by  altering  the  setting  on  the  plenum  valves, 
the  plenum  valve  settings  could  be  held  constant.  Only  the 
stop  valve  was  needed  to  initiate  and  terminate  a  particular 
experimental  run. 

In  anticipation  of  testing  in  the  Naval  Postgraduate 
School  smoke  tunnel,  a  two- inch  pipe  was  channeled  to  the 
vicinity  of  the  smoke  tunnel  from  the  storage  tanks.   In 
parallel  with  the  outlet  in  the  static  testing  area,  the 
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pipe  routed  to  the  smoke  tunnel  was  sealed  during  static 
testing.   Simple  transition  to  smoke  tunnel  testing  can  be 
performed  by  sealing  the  outlet  in  the  static  testing  area 
and  moving  the  regulator  and  stop  valves  to  the  smoke  tun- 
nel outlet. 

The  plenum  has  eight  outlets  of  which  only  six  are  needed 
to  supply  the  ejectors  on  the  wing  and  canard.   Three-quarter- 
inch-diameter  automotive  heater  hoses  were  used  to  connect 
the  plenum  to  the  Coanda  slots  for  both  wing  and  canard. 
One-inch  hoses  were  needed  for  both  center  ejectors  (Fig. 
5). 

B.   INLET  SIMULATION 

In  order  to  depict  adequately  aircraft  operation  in  the 
static  position,  some  means  of  duplicating  the  effect  of 
the  inlets  needed  to  be  accomplished.   Upon  commencing  work 
on  the  XFV-12A,  only  the  main  inlet  had  been  fabricated 
with  no  means  of  providing  suction  for  inlet  flow. 

The  auxiliary  inlet  was  fashioned  from  copper  sheet  metal 
and  attached  to  the  top  of  the  aircraft  model  as  shown  in 
Figure  6.   For  both  the  main  and  auxiliary  inlets,  ducts 
were  constructed  for  the  purpose  of  establishing  the  correct 
inlet  flow  rate.   The  ducts  are  shown  in  Figure  7.   A  suf- 
ficient quantity  of  inlet  flow  is  pumped  to  satisfy  the 
continuity  equation  as  applied  to  inlet  and  exit  conditions. 
Continuity  required  that  the  sum  of  the  mass  flow  rates 
entering  the  main  and  auxiliary  inlets  was  equal  to  the  mass 
flow  rate  of  primary  air  from  the  wing  and  canard.   (For 
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purposes  of  these  experiments,  the  contribution  of  fuel  to 
the  mass  flow  rates  of  the  exhaust  jets  was  assumed  to  be 
negligible . ) 

A  first  attempt  at  entraining  ambient  air  through  the 
inlets  was  done  using  the  ejector  principle.   Copper  tubing 
with  six  ejectors  was  fabricated  to  blow  air  along  the  main 
inlet  channel  near  the  face  of  the  inlet.   A  photograph  of 
the  ejector  appears  in  Figure  9.   A  similar  copper  ejector 
was  constructed  for  the  same  purpose  in  the  auxiliary  inlet. 
These  ejectors  were  found  to  function  but  could  not  provide 
sufficient  mass  flow  to  match  that  leaving  the  Coanda  slots 
and  ejectors  of  the  propulsive  system.   Thus,  some  other 
means  of  duplicating  the  suction  effect  at  the  inlets  had 
to  be  devised. 

A  suction  blower  was  attached  to  the  inlet  channels  and 
sealed  so  that  all  of  the  entrained  air  was  expelled  out  the 
exit  plane  of  the  blower;  see  Figure  10.   Using  this  method 
the  dynamic  pressure  measured  at  the  exit  plane  of  the 
blower  was  found  to  be  0.6  inches  of  water.   For  purposes  of 
analysis  of  the  continuity  equation,  the  flow  was  assumed  to 
be  incompressible  between  the  exit  plane  of  the  blower  and 
the  inlets.   With  this  simplification,  velocities  and  cor- 
responding dynamic  pressures  at  the  inlet  and  the  blower 
exit  plane  were  dependent  only  on  the  relative  areas. 

U,A,  =  U  A  (9) 

4  4mm  v  ' 

and 

q-     U,        A 

qm      m        A4 
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The  dynamic  pressure  at  the  exit  plane  of  the  motor  dic- 
tates the  conditions  which  must  exist  at  the  ejectors  for 
continuity  to  be  satisfied  for  the  system. 

Upon  comparison  of  geometries,  the  required  velocity  at 
the  main  inlet  was  found  to  be  27.5  ft/sec.   The  correct 
velocity  at  the  ejectors  and  Coanda  slots  of  the  wing  and 
canard  was  then  136  ft/sec.   Applying  the  augmentation 
ratio  determined  previously,  the  mean  velocity  at  the  exit 
plane  of  the  diffusers  must  be  26.8  ft/sec. 

C.   MODEL  MOUNTING 

Once  the  inlets  had  been  constructed,  the  model  was 
mounted  on  three-quarter-inch  plywood  so  that  a  plane  of 
symmetry  was  established  for  testing  purposes.   The  plane 
of  symmetry  is  commonly  known  as  a  reflection  plane.   The 
other  side  of  the  reflection  plane  was  cut  away  at  the  in- 
terior of  the  model  to  allow  for  the  introduction  of  air 
from  the  plenum.   A  photograph  of  the  reflection  plane 
mounting  at  the  static  testing  site  is  shown  in  Figure  8. 

The  ground  was  simulated  by  another  piece  of  plywood 
mounted  perpendicular  to  the  model  plane  of  symmetry  and 
parallel  to  the  ground.   This  piece  was  attached  to  the 
reflection  plane  with  bolts  and  could  be  moved  up  and  down 
with  respect  to  the  aircraft  at  one-inch  intervals.   The 
upper  limit  was  analogous  to  placing  the  aircraft  on  the 
ground  with  the  gear  lowered.   The  lower  limit  simulated 
the  aircraft  at  an  altitude  of  20  feet  and  was  accomplished 
by  removing  the  ground  plane  entirely. 
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As  the  "ground'1  was  moved  down  from  the  aircraft  during 
testing,  it  became  necessary  to  use  a  larger  piece  of  ply- 
wood to  compensate  for  the  enlarged  area  which  was  under 
the  influence  of  the  aircraft's  exhaust  jets.   This  was  done 
by  placing  a  larger  piece  of  plywood  covering  52  square  feet 
directly  over  the  original  "ground." 
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V.   PRELIMINARY  EXPERIMENTAL  WORK  AND  ANALYSIS 

A.   NATURE  OF  RECIRCULATION 

Reingestion  of  hot  exhaust  gases  is  a  problem  of  particu- 
lar interest  to  designers  of  VTOL  aircraft  and  STOL  aircraft 
when  operating  at  forward  speeds  of  less  than  25  knots  [Ref. 
12].   This  phenomenon  involves  the  entrainment  of  some  por- 
tion of  one  or  more  exhaust  jets  into  the  intakes  of  the 
lift  engines.   Adverse  effects  on  performance  include  a  de- 
crease in  the  available  thrust,  large  increase  in  inlet 
temperature,  and  distortion  in  the  temperature  profile  at 
the  inlets  [Ref.  7] .   These  last  two  effects  can  lead  to 
stalling  of  the  compressor  and  engine  surge  [Ref.  7] .   Ver- 
tical takeoff  and  landing  operations  over  areas  with  loose 
terrain  also  pose  a  problem  of  ingestion  of  particles  kicked 
up  by  the  exhaust  jets  [Ref.  13] .   Potential  loss  of  avail- 
able control  power  was  illustrated  dramatically  by  a  XC-142. 
A  landing  accident  in  1965  was  attributed  to  a  loss  in  con- 
trol power  due  to  recirculation  of  hot  exhaust  gases  [Ref. 
12]  . 

Exhaust  gas  reingestion  has  been  classified  according  to 
the  method  in  which  exhaust  gas  is  transported  to  the  vicinity 
of  the  inlet.   The  first  type  is  a  diffusion  of  the  exhaust 
gases  into  the  ambient  air.   The  buoyant  force  of  the  low 
density  exhaust  carries  it  up  from  the  ground  to  the  vicinity 
of  the  inlet.   This  is  commonly  called  the  "buoyant"  effect 
[Ref.  5]. 
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The  second  type  of  recirculation  involves  an  inter- 
action of  the  exhaust  jets  with  the  ground  and  subsequently 
with  one  another  so  that  exhaust  gases  are  forced  upward  in 
a  "fountain"  effect  between  the  two  exhaust  jets.   This  type 
usually  requires  less  time  for  its  effects  to  be  felt.   The 
fountain  effect  usually  is  considered  to  be  the  most  serious 
of  the  two  types  in  terms  of  potential  hazard  to  aircraft 
operation  [Ref.  5]. 

The  reingestion  problem  is  one  which  cannot  be  treated 
analytically  at  present.  Rather,  experimental  results  are 
the  only  reliable  source  of  information  concerning  the  re- 
circulation characteristics  of  a  particular  airframe  shape 
and  engine  configuration.  Even  the  interpretation  of  re- 
sults has  resulted  in  conflicting  conclusions.  For  example, 
independent  tests  arrived  at  opposite  conclusions  with  re- 
gard to  the  most  significant  type  of  reingestion  transport 
[Refs.  5,7]. 

B.   ANALYSIS  OF  VALIDITY  OF  SMALL-SCALE  TESTING 

In  order  for  experimental  results  to  be  valid,  reinges- 
tion characteristics  observed  for  the  scale  model  must  be 
applicable  to  the  full-scale  aircraft.   Also,  since  operating 
conditions  of  the  aircraft  may  differ  from  those  imposed  on 
the  model  during  testing,  reingestion  must  be  evaluated  as 
a  function  of  differing  exit  temperature  and  pressure 
conditions . 

In  an  investigation  of  scaling  effects  on  VTOL  recircu- 
lation, recirculation  was  determined  to  be  scalable  both 
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statically  and  dynamically  [Ref.  6].   (The  experimental 
variables  used  were  temperature,  temperature  gradient,  and 
inlet  distortion.)   Another  analysis  determined  reingestion 
to  be  independent  of  the  temperature  and  pressure  profiles 
exhibited  by  the  engine  over  the  range  tested  [Ref.  4] . 
These  conclusions  indicated  that  the  experimental  test  con- 
ditions provided  significant  results  even  though  full  scale 
operating  conditions  were  not  duplicated. 

With  specific  regard  to  the  XFV-12A,  an  evaluation  of 
the  characteristics  of  thrust  augmentation  indicated  that 
the  augmentation  ratio  of  a  system  is  almost  independent  of 
nozzle  pressure  ratio  [Ref.  14].   Consequently,  conditions 
unique  to  propulsive  lift  employing  thrust  augmentation  were 
duplicated  at  the  test  conditions. 

C.   SMOKE  AND  LASER  TESTING 

The  first  experiments  performed  to  determine  the  existence 
of  recirculating  exhaust  gases  were  done  using  smoke  intro- 
duced into  the  flow  and  a  portable  laser  for  detection.   A 
smoke  jet  was  produced  from  an  oil  smoke  generator  similar 
to  those  used  in  the  NPS  smoke  tunnel.   Using  an  air  flow 
from  one  of  the  extra  outlets  of  the  operating  plenum,  the 
smoke  was  forced  through  a  narrow  length  of  plastic  tubing 
to  the  vicinity  of  the  canard  ejector. 

The  smoke  was  laminar  for  about  2  or  3  inches  upon 
leaving  the  tube,  followed  by  transition  to  turbulence  as 
it  was  entrained  with  secondary  air  into  the  canard.   Sub- 
sequent to  entrainment  the  smoke  was  impossible  to  trace 
without  the  use  of  the  laser. 
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The  smoke  experiments  were  successful  in  demonstrating 
a  small  degree  of  recirculation,  although  no  quantitative 
results  could  be  determined.   They  also  demonstrated  that 
the  main  and  auxiliary  inlet  ejectors  functioned  as  desired 
but  were  unable  to  develop  sufficient,  amounts  of  inlet  flow 
to  match  the  wing  and  canard  ejectors.   The  results  of  smoke 
testing  were  inconclusive;  however,  the  results  did  indicate 
that  further  investigation  was  desirable. 

D.   EXPLORATORY  FREON  TESTING 

Experiments  conducted  in  similar  fashion  to  the  smoke 
tests  were  performed  using  an  aerosol  can  containing  freon 
instead  of  smoke.   A  freon  leak  detector  normally  used  to 
detect  refrigerant  leaks  was  used  to  measure  reingestion. 

These  results  confirmed  that  exhaust  gases  were  being 
re  ingested,  though  there  was  again  inadequate  suction  at 
both  intakes.   This  experiment  demonstrated  the  need  for  a 
better  method  of  inlet  simulation  as  well  as  further  analysis 
of  reingestion  characteristics. 
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VI.   REINGESTION  EXPERIMENTS 

A.   DESCRIPTION  OF  LEAK  DETECTOR 

As  a  detection  device,  a  General  Electric  refrigerant 
leak  detector  was  acquired.   Normally  used  to  detect  leaks 
as  small  as  1/2  oz/yr  in  coolant  systems,  it  was  aptly 
suited  for  use  in  detecting  exhaust  gases  contaminated  with 
freon.   The  detector  is  illustrated  in  Figure  11. 

The  sensitivity  of  the  detector  was  varied  by  a  knob 
on  the  face  of  the  instrument  and  by  a  switch  which  indi- 
cated either  a  high  or  low  sensitivity  setting.   The  probe 
held  a  light  which  flashed  when  freon  was  absorbed.   The 
flash  rate  gave  a  qualitative  indication  of  the  absorption 
rate,  although  the  flash  frequency  was  not  proportional  to 
the  freon  contamination  of  the  mass  flow  being  absorbed. 

A  moderate  rate  of  reingestion  was  defined  as  that 
flash  frequency  which  was  equal  to  that  obtained  when  the 
probe  was  placed  directly  over  the  reference  leak.   Since 
the  reference  leak  has  a  known  mass  flow  rate,  it  can  be 
used  as  a  reference  value.   Three  other  ordinal  quantitative 
descriptions,  heavy,  light,  and  trace,  were  used  to  repre- 
sent deviations  from  the  reference  "moderate"  condition. 

Before  each  experimental  run,  the  detector  was  adjusted 
by  reducing  or  increasing  the  sensitivity  until  at  free- 
stream  conditions,  the  probe  light  just  ceased  to  flash. 
Thus  each  set  of  data  was  acquired  using  the  same  reinges- 
tion descriptions.   Contamination  of  the  ambient  air  at  low 
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levels  of  freon  forced  experiments  to  be  performed  at  sub- 
stantial time  intervals.   This  allowed  the  background 
freon  contamination  to  dissipate  and  provided  time  to  re- 
supply  the  compressed  air  storage  tanks. 

B.   EXPERIMENTAL  DEPENDENT  VARIABLES 

Of  interest  in  performing  the  experiments  was  the  flow 
path  of  recirculating  gases  from  the  exit  plane  of  the  dif- 
fuser  of  both  the  canard  and  wing  to  the  main  and  auxiliary 
inlets.   It  was  anticipated  that  these  recirculation  char- 
acteristics would  vary  with  height  above  the  ground  and 
with  the  angle  between  the  chord  of  the  wing  and  canard 
flaps  and  the  ground.   The  investigation  covered  recircula- 
tion characteristics  as  height  above  the  ground  varied  from 
zero  to  about  24  inches  (20  feet  in  full-scale  dimensions). 
The  angle  of  the  canard  and  wing  flaps  was  varied  for  each 
height  at  10  degree  increments  from  the  vertical  to  an 
angle  of  30  degrees  from  the  vertical.   Thus  the  pattern  of 
reingestion  was  investigated  for  the  aircraft  operating  in 
both  the  STOL  and  VTOL  modes. 

In  addition,  a  profile  of  the  distribution  of  recircula- 
ted exhaust  at  the  main  inlet  was  determined.   This  was  done 
by  converting  the  six  main  inlet  ejectros  previously  employed 
unsuccessfully  to  the  role  of  gas  sampling  ports.   Tests  were 
performed  at  a  height  above  ground  of  zero  and  with  the  wing 
and  canard  flaps  in  the  VTOL  configuration  only.   Profiles 
were  acquired  which  reflected  the  distribution  of  exhaust 
gases  from  the  front  or  rear  Coanda  slots  or  the  ejector  in 

the  canard. 
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C.   QUANTITATIVE  RE INGEST ION  ANALYSIS 

Some  means  were  necessary  of  quantitatively  specifying 
the  meaning  of  the  descriptions  of  degree  of  reingestion. 
Since  the  term  "moderate"  has  been  assigned  to  that  amount 
of  intake  corresponding  to  the  reference  leak,  it  provided 
the  only  basis  for  a  quantitative  analysis. 

A  result  can  be  obtained  analytically  in  a  straightfor- 
ward manner  if  some  assumptions  are  made  concerning  the  con- 
ditions under  which  the  experiments  were  conducted  and  the 
nature  and  properties  of  the  freon.   First,  it  was  assumed 
that  thermodynamic  properties  could  be  evaluated  at  70°F. 
The  freon  was  assumed  to  have  the  molecular  weight  of  Freon- 12 
(CC1~F~) ;  see  [Ref .  15] .   The  freon  in  the  aerosol  can  was 
assumed  to  be  a  mixture  of  vapor  and  liquid  at  saturation 
conditions . 

Five  mass  flow  rates  must  be  calculated  using  known 
properties  of  freon  and  air  and  the  previously  outlined  as- 
sumptions.  The  mass  flow  rate  of  freon  from  the  reference 
leak  is  required;  this  is  equal  to  the  mass  flow  rate  for 
moderate  recirculation.   The  mass  flow  rate  of  air  in  the 
diffuser  was  also  needed,  as  well  as  the  mass  flow  rate  of 
the  detection  probe.   Finally,  the  mass  flow  rates  of  the 
inlets  and  the  aerosol  can  were  calculated.   The  various 
mass  flow  rates  and  corresponding  sources  and  sinks  are 
diagramatically  shown  in  Figure  12. 

The  mass  flow  rate  of  the  reference  leak  is  known  to  be 
1/2  oz/yr  from  the  instructions  for  use  of  the  instrument. 
When  the  units  are  altered  to  yield  lbm/sec,  the  result  is 
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m  =  0.99  x  10"9  lbm/sec. 
o 

The  mass  flow  rate  of  the  aerosol  can  was  resolved  using 

the  momentum  conservation  equation  and  solving  for  velocity. 

u| .  am  .  (id 

The  pressures  and  density  were  obtained  from  a  table  of 
properties  of  saturated  freon  [Ref.  15].   The  calculated 
velocity  is 

U1  =  567  ft/sec. 

The  diameter  of  the  tube  leading  from  the  aerosol  can  was 
measured  to  be  0.025  inches.  The  mass  flow  rate  was  then 
readily  calculated. 

m1  =  3.92  x  10    lbm/sec. 

The  mass  flow  rate  in  the  diffuser  was  computed  using 
the  results  obtained  for  the  augmentation  ratio  [Ref.  8] 
and  the  exit  conditions  of  the  wing  and  canard  diffusers. 
Using  a  primary  air  velocity  of  136  ft/sec,  the  mass  flow 
rate  of  the  canard  diffuser  is 

m2  =  4.58  x  10" 1    lbm/sec. 

The  mass  flow  rate  in  the  probe  required  a  supplementary 
experiment.   The  probe  was  attached  to  a  manometer  at  one 
end.   The  other  end  led  to  a  reservoir  of  fluid  whose  upper 
surface  was  exposed  to  ambient  conditions.   For  small  changes 
in  pressure,  incompressibility  can  be  assumed  for  the  air. 
The  flow  rate  of  gas  being  drawn  into  the  probe  could  be  in- 
ferred by  observing  the  rate  of  displacement  of  the  manometer 

fluid. 
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Using  the  speed  of  the  air  in  the  manometer  tube,  the 
mass  flow  rate  of  the  gas  in  the  detector  probe  could  be 
computed  directly  from  the  continuity  equation  for  incom- 
pressible flow.   Subsequent  applications  of  these  results 
assumed  that  air  was  the  dominant  gas  being  absorbed  by  the 
probe . 

U  A   =  uYA..  (12) 

p  p    M  M  v  J 

Solving  for  the  only  unknown,  U  ,  yields  a  velocity  of  gas 

in  the  probe  of  5.37  in/sec.   The  mass  flow  rate  of  the 

probe  was 

m_  =  4.82  x  10    lbm/sec. 

The  mass  flow  rate  of  the  main  inlet  was  calculated 
using  measurements  taken  at  the  exit  plane  of  the  electric 
blower.   The  dynamic  pressure  was  recorded  which  yielded 
the  velocity  at  the  blower  exit.   Comparing  geometries  from 
the  analysis  described  in  Section  IVB,  the  properties  at 
both  inlets  could  be  calculated.   The  resultant  mass  flow 
rate  for  the  main  inlet  was 

m.    =  6.29  x  10    lbm/sec. 

The  mass  flow  rate  of  the  auxiliary  inlet  was  assumed 
to  be  the  same  as  that  of  the  main  inlet.   The  area  of  the 
auxiliary  inlet  is  greater  than  that  of  the  main  inlet. 
However,  two  ducts  of  equal  area  connected  the  inlets  to 
the  blower.   These  ducts  were  smaller  in  cross  sectional 
than  either  inlet  and  consequently  dictated  the  mass  flow 
rate  of  each.  »  " 
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The  first  problem  in  quantitative  analysis  was  the  com- 
putation of  the  fraction  of  mass  flow  in  the  detection  probe 
that  was  freon.   For  a  moderate  reingestion  rate,  this  was 
represented  by  the  ratio  of  the  mass  flow  rate  in  the  refer- 
ence leak  to  the  mass  flow  rate  of  the  probe. 

F1  =  m  /m3  =  2.06  x  lO-*. 

The  second  problem  was  finding  the  concentration  of 
freon  in  the  diffuser  of  the  canard.  Since  this  concentra- 
tion was  assumed  to  be  constant  from  the  point  of  introduc- 
tion to  the  ejector,  it  was  found  by  computing  the  ratio  of 
mass  flow  rate  of  freon  introduced  into  the  air  hose  to  the 
mean  flow  rate  of  air  at  the  exit  plane  of  the  diffuser. 

F2  =  m1/in2  =  8.58  x  10-3. 

A  reingestion  percentage  was  then  calculated  yielding 
the  concentration  of  freon  being  ingested  as  compared  to 
that  being  expelled.   For  moderate  reingestion  conditions 

F  =  F1/F2  =  2.40%. 

Consequently  a  moderate  reading  indicated  that  2.40%  of  the 
diffuser  exit  mass  flow  was  being  reingested  by  the  main 
inlet . 

Further,  it  was  desirable  to  know  the  percentage  of 
inlet  flow  that  was  recirculated  from  the  engine  exhausts. 
This  was  computed  simply  by  multiplying  the  reingestion 
percentage  by  the  mass  flow  rate  of  the  diffuser,  then 
dividing  by  the  mass  flow  rate  of  the  main  inlet. 

G  =  F  x  m2/m4  =  17.8%. 
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Therefore  a  moderate  reingestion  reading  indicated 
that  about  18%  of  the  inlet  flow  was  comprised  of  recircu- 
lated exhaust  gases.   Here  exhaust  gases  were  defined  as 
the  mixture  of  entrained  secondary  air  and  primary  air  ex- 
pelled from  the  canard  or  wing  augmentor. 

D.   EXPERIMENTAL  PROCEDURE 

Freon  was  injected  through  drilled  holes  in  the  heater 
hoses  leading  to  the  Coanda  slots  and  the  ejectors.   The 
aerosol  can  could  be  moved  readily  to  introduce  freon  into 
any  given  flap.   The  six  points  of  introduction  were  as- 
signed station  numbers  sequentially  from  the  forwardmost 
canard  flap  aft.   The  flaps  and  their  accompanying  station 
numbers  are  shown  in  Figure  13.   Holes  were  drilled  in  the 
ducts  leading  from  the  inlets  to  which  the  detector  probe 
was  attached.   The  wing  and  canard  flaps  were  set  by  mea- 
suring the  angle  between  the  chord  of  each  flap  and  the 
vertical.   The  stiffness  of  the  hoses  was  sufficient  to 
maintain  the  position  of  the  wing  and  canard  flaps  through 
out  a  set  of  experiments. 

The  experiments  were  conducted  by  varying  the  height 
above  the  ground  with  a  constant  wing  and  canard  angle. 
Measurements  were  taken  at  both  the  main  and  auxiliary  in- 
lets for  each  of  the  six  points  of  introduction.   Thus 
twelve  readings  were  taken  at  each  height  above  the  ground 
for  each  wing  and  canard  angle.   The  wing  and  canard  were 
rotated  together,  so  that  there  were  no  experimental  runs 
for  which  the  wing  and  canard  angles  were  different.   A 
total  of  288  reingestion  readings  were  taken. 
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The  results  were  plotted  according  to  the  twelve  com- 
binations of  inlet  and  point  of  introduction.   The  twelve 
reingestion  profiles  are  shown  in  Figures  14-25.   The 
abscissa  of  each  graph  is  the  angle  between  the  chord  of  the 
wing  and  canard  flaps  and  the  ground.   The  ordinate  is  the 
height  above  the  ground,  non-dimensionalized  with  respect 
to  a  canard  reference  length.   The  canard  reference  length 
was  defined  as  the  square  root  of  the  area  of  the  entrance 
plane  of  the  canard  augmentor.   Its  value  is  4.90  in. 

The  experiments  covering  the  inlet  reingestion  profile 
were  performed  with  the  canard  since  the  wing  did  not  con- 
tribute substantially  to  main  inlet  recirculation.   Five  of 
the  six  different  gas  sampling  ports  in  the  inlet  face  were 
sealed  with  adhesive  tape  at  any  one  time,  so  that  the  rein- 
gestion reading  measured  only  streamtubes  passing  near  the 
open  port.   The  response  of  all  six  sampling  ports  on  the 
inlet  face  was  recorded  for  each  of  the  three  points  of  in- 
troduction.  A  total  of  18  readings  were  taken.   Patterns  of 
the  main  inlet  reingestion  are  shown  in  Figures  26-28. 
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VII.   EXPERIMENTAL  RESULTS  AND  INTERPRETATION 

A.   POTENTIAL  FLOW  SOLUTION 

The  digital  potential  flow  solution  yielded  a  result 
which  indicated  that  "supercirculation"  was  present  at  the 
rearmost  canard  flap.   Here  supercirculation  was  defined  to 
occur  when  a  streamline  circumscribed  any  particular  lifting 
surface.   Supercirculation  was  present  throughout  the  range 
of  variables  of  velocity  in  the  diffuser  and  mass  flow  dis- 
tribution entering  the  canard  augmentor.   There  was  also  a 
rapid  return  of  the  exhaust  profile  to  freestream  conditions 
only  a  few  chord  lengths  aft  of  the  canard. 

The  electrostatic  analogy  solution  eliminated  inaccuracies 
in  the  coarse  mesh  size.   The  flow  around  the  canard  was  ab- 
sent of  supercirculation.   However,  a  stagnation  point  was 
produced  on  the  upper  surface  of  the  rear  canard  flap.   The 
influence  of  the  wing  and  its  exhaust  was  apparent  in  its 
influence  on  the  canard  flow.   The  tendency  of  the  canard 
exhaust  jet  to  return  to  freestream  conditions  was  negated 
by  the  presence  of  the  wing  and  its  flow  field.   The  highest 
velocities  in  the  flow  field  occurred  between  the  wing  and 
the  simulated  ground.   The  bending  of  the  exhaust  jets  and 
resultant  vortex  shedding  was  evident  primarily  around  the 
wing.   Supercirculation  was  evident  around  the  rearmost  flap 
of  the  wing  for  this  case. 

Using  the  results  of  the  electrostatic  analogy  as  an 
indication  of  the  velocity  and  shape  of  the  exhaust  jets, 
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an  analysis  was  done  to  determine  the  effect  of  the  exhaust 
jets  and  vortex  shedding  on  the  lift  of  the  aircraft.   Using 
the  results  shown  in  Figure  4,  the  radii  of  curvature  of  the 
exhaust  streamlines  and  the  distances  between  the  adjacent 
streamlines  in  the  normal  direction  were  measured  graphically. 
This  was  done  in  the  vicinity  of  the  trailing  edge  of  the  rear 
wing  flap. 

The  rate  of  change  of  pressure  in  the  normal  direction 
is  related  to  the  local  dynamic  pressure  and  streamline  radius 
of  curvature  by 

U2      8P  ri,. 

P    n =   -  ^—  •  (13) 

Ka  R. .      dn  v  J 

Considering  streamlines  one  and  two  of  Figure  4,  the  above 
equation  can  be  integrated  to  yield  the  following  result. 

P2  "  Pl  "  "  Pa  R^  ("2  -  V-  W 

Nondimensionalizing  the  pressure  difference  with  respect  to 
the  freestream  dynamic  pressure,  the  following  relation  can 
be  derived. 

(P2  -  P  )      2(AnJ2 

—  =  "  d 7 T  •  (15) 

qco  Ri2(n2-ni} 

The  results  of  these  measurements  and  a  resultant  pressure 
coefficient  across  the  exhaust  jet  are  shown  in  Appendix  C. 
Further,  since  the  lift  coefficient  can  be  represented 
by  the  difference  in  the  average  pressure  at  the  top  of  a 
surface  and  the  bottom  of  the  surface  divided  by  the  dynamic 
pressure,  the  resultant  value  for  AP/q^   across  the  jet  is 
analogous  to  a  lift  coefficient.   The   numerical  value  for 
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the  pressure  coefficient  across  the  jet  is  greater  than  that 
which  would  have  been  obtained  from  a  three-dimensional  ap- 
plication of  potential  flow. 

B.   RECIRCULATION 

Recirculation  was  most  severe  in  the  VTOL  configuration 
when  the  aircraft  was  situated  on  the  ground.   The  main  in- 
let was  the  largest  recipient  of  exhaust  gases.   In  general, 
reingestion  decreased  as  height  above  the  ground  increased. 
However,  for  some  cases,  reingestion  increased  as  height 
aobve  the  ground  increased  and  then  decreased  again  as  the 
height  was  further  increased;  for  example,  see  Figure  15, 
9  =  60°,  and  Figure  16 ,  G  =  90°.   The  highest  readings  re- 
sulted when  the  point  of  introduction  was  at  the  forward 
Coanda  slot  of  the  canard. 

The  results  of  introducing  freon  into  the  canard  or 
wing  were  not  independent  of  the  point  of  introduction. 
Usually  there  was  considerable  difference  in  the  response 
from  station  to  station  within  the  wing  or  canard.  For 
example,  see  Figures  14,  15,  and  16  for  the  response  of 
the  main  inlet  to  introduction  points  within  the  canard. 

The  angle  that  the  wing  and  canard  made  with  the  ground 
had  considerable  influence  on  recirculation.   With  the  wing 
and  canard  at  an  angle  of  60  degrees  with  the  ground,  the 
wing  was  never  a  source  of  reingestion.   The  canard  pro- 
vided some  reingested  gases  at  all  altitudes  and  configura- 
tions.  Deviations  in  wing  and  canard  angle  with  the 
vertical  resulted  in  a  substantial  reduction  of  reingestion 
for  all  altitudes  and  configurations. 
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The  distribution  tests  at  the  main  inlet  indicated  that 
reingestion  was  concentrated  most  heavily  on  the  side  most 
distant  from  the  plane  of  symmetry  of  the  aircraft  when  the 
point  of  introduction  was  station  1;  see  Figure  13  for 
station  identification  and  Figure  26  for  the  measured  re- 
sult.  For  station  2,  however,  the  reingestion  properties 
were  most  pronounced  at  the  upper  portion  of  the  main  in- 
let; see  Figure  27.   When  freon  was  introduced  into  station 
3,  there  was  an  even  distribution  over  the  entire  main  inlet 
face;  see  Figure  28.   The  level  of  reingestion  was  only 
trace  for  this  case. 

The  response  time  for  all  reingestion  readings  was 
small,  indicating  that  the  fountain  effect  was  the  only  ef- 
fect contributing  to  recirculation.   There  was  no  opportunity 
for  the  "buoyant"  effect  to  be  felt,  however,  since  there 
was  no  significant  temperature  difference  between  the  pri- 
mary and  ambient  air.   Responses  occurring  a  considerable 
length  of  time  after  the  introduction  of  freon  were  attri- 
buted to  diffusion  and  mixing  of  freon  into  the  ambient 
air.   These  responses  were  never  greater  in  magnitude  than 
trace. 
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VIII.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  potential  flow  solution  was  valuable  in  obtaining  an 
approximation  to  the  pressure  coefficient  across  the  exhaust 
jet.   This  coefficient  provided  insight  into  the  contribu- 
tion of  a  two-dimensional  propulsive  lift  jet  to  the  flow 
field  of  the  aircraft. 

Reingestion  poses  a  problem  primarily  to  the  main  inlet 
under  VTOL  conditions.   The  magnitude  of  the  problem  was 
greater  than  the  181  obtained  with  moderate  reingestion 
for  all  altitudes.   The  canard  was  the  most  significant 
source  of  reingested  exhaust  gases,  and  was  the  only  source 
to  supply  readings  above  moderate.   The  wing  contributed  no 
gases  for  reingestion  in  configurations  other  than  VTOL. 

B.  RECOMMENDATIONS 

Further  work  would  be  desirable  to  determine  the  effect 
of  a  forward  velocity  as  simulated  in  the  wind  tunnel  on 
reingestion  properties.   Closer  scrutiny  might  be  applied 
to  the  flow  path  from  the  canard  ejector  and  Coanda  slots 
to  the  main  inlet. 
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APPENDIX  A:   SAMPLE  DERIVATION  OF  DIFFERENCE  EQUATION 

Each  linear,  algebraic  equation  was  derived  as  follows 

iKx,y+k) 

5}  k 

*(x-j,y) '  ^'        '^(x+j ,y) 

^(x,y-k) 

beginning  with  Taylor  series  for  each  axis  of  the  domain. 

♦(x.y+k)  =  *(x,y)  +  ^y   (x,y)  +  j£  0  (x,y)  +  ... 
and 

i  r  ^^     ,   r  ^     k  3  ^   r     ^_,_k282^,     -> 

iKx,y-k)  =  *(x,y)  -  ^-  (x,y)  +  JT  df*    (  ,y^  "  •" 

Adding  the  above  equations  and  solving  for  the  second  par- 
tial, the  following  equation  results. 

f^  =  |r  [  <Hx,y+k)  -  2*(x,y)  +  *(xfy-k)  ]  +  0(4). 

Similar  analysis  for  the  second  partial  in  the  x-direction 
yields 

|1|  =  |_  f  ^(x+j>y)  .  2^(x,y)  +  if^(x-j,y)  ]  +  0(4). 

Applying  Laplace's  equation 

ili  +  *1±   =  n 

8X2      3y2 

each  solution  point  is  represented  as  a  linear  function  of 
its  four  adjacent  points. 

i2  k2 

<Kx,y)  =  2(j2+T^J  t  <Kx>y+k)  +  ip(x,y-k)]  +  JJJ^W) 

[  ^(x+j  ,y)  +  i|i(x~j  ,y)]  . 
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APPENDIX  B:   MASS  FLOW  RATE  CALCULATIONS  OF  QUANTITATIVE 
ANALYSIS 


The  following  data  was  obtained  either  experimentally 

or  from  tables  of  known  properties  of  freon  or  air. 

p   =  7.66  x  10"2  lbm/ft3 
a 


pf 

= 

2.028 

lbm/ft3 

ul 

= 

567 

ft/sec 

U2 

= 

136 

ft/sec 

U 
P 

= 

0.45 

ft/sec 

U4 

= 

27.5 

ft/sec 

Al 

= 

4.91  x 

lO"" 

in2 

A2 

= 

0.807 

2 

A 
P 

= 

.0204 

in2 

A4 

= 

4.3 

in2 

The  mass  flow  rates  used  in  the  quantitative  analysis 
were  calculated  using  the  above  data  and  appropriate  conver 
sion  factors  to  maintain  consistency  in  units. 


ml  =  P£U!A!  =  3-92  x  10  lbm/sec 

m2  _  paU2A2  =  4'58  x  10_1  lbm/sec 

m-  =  P  U  A   =  4.82  x  10"6  lbm/sec 
j          a  p  p 

m4  =  PU.A4  =  6.29  x  10-2  lbm/sec 
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APPENDIX  C:   RESULTS  OF  POTENTIAL  FLOW  LIFT  JET  ANALYSIS 


Streamtube 


An 


n .  -n. 


R 


2(An  2) 


ii     R. . fn.-n.] 


1,2 


.39 


.25 


.75 


4.16 


2,3 


.19 


.125 


1.20 


2.53 


3,4 


.19 


.156 


1.69 


-  1.44 


4,5 


.19 


.31 


2.38 


0.53 


2  (An/) 

R.  .  (n.-n.) 
il  3      i 


=  8.66 


Note:   Velocity  of  flow  at  augmentor  exit  planes  was  twice 
that  of  freestream  conditions. 
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(a)   ORDINARY  FINITE  SPAN  WING 


(b)   JET  FLAP  FINITE  SPAN  WING 


Figure  1.   SHED  VORTICITY  ASSOCIATED  WITH  LIFTING  WINGS 
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Note:   Resultant  stream  function  values  are  shown  in  sample 
computer  program  output. 


Figure  2.   DISCRETIZED  DOMAIN 
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Figure  4.   ELECTROSTATIC  ANALOGY  RESULTS 
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ENTRANCE  PLANE 


TO 
AMBIENT 
AIR 


CANARD  DIFFUSER 
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SECONDARY 
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FROM  AMBIENT 
AIR 


PRIMARY  AIR 


FROM 
AEROSOL 
CAN 


Figure  12.   PATH  OF  REINGESTED  EXHAUST  GASES. 
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(a)   CANARD  STATIONS 


\ 

5 


(b)   WING  STATIONS 
Figure  13.   IDENTIFICATION  OF  WING  AND  CANARD  FLAPS 
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Figure  14.   REINGESTION  PATTERN;  MAIN  INLET,  STATION  1 
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Figure  15.   REINGESTION  PATTERN;  MAIN  INLET,  STATION  2 
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Figure  16.   REINGESTION  PATTERN;  MAIN  INLET,  STATION  3 
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Figure  17.   REINGESTION  PATTERN;  MAIN  INLET,  STATION  4 
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Figure  18.   REINGESTION  PATTERN;  MAIN  INLET,  STATION  5 
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Figure  19.   REINGESTION  PATTERN;  MAIN  INLET,  STATION  6 
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Figure  20.   REINGESTION  PATTERN;  AUXILIARY  INLET,  STATION  1 
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Figure  21.   REINGESTION  PATTERN;  AUXILIARY  INLET,  STATION  2 
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Figure  22.   REINGESTION  PATTERN;  AUXILIARY  INLET,  STATION  3 
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Figure  23.   REINGESTION  PATTERN;  AUXILIARY  INLET,  STATION  4 
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Figure  24.   REINGESTION  PATTERN;  AUXILIARY  INLET,  STATION  5 
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Figure  25.   REINGESTION  PATTERN;  AUXILIARY  INLET,  STATION  6 
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O  SAMPLE  PORT 


Figure  26.   REINGESTION  PATTERN;  MAIN  INLET  FACE,  STATION  1 
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O  SAMPLE  PORT 


Figure  27.   REINGESTION  PATTERN;  MAIN  INLET  PACE,  STATION  2 
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O  SAMPLE  PORT 


Figure  28.   REINGESTION  PATTERN;  MAIN  INLET  FACE,  STATION  3 
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COMPUTER  PROGRAM 

On  the  next  page,  the  text  of  the  computer  program 
begins  which  solved  Laplace's  equation  using  finite  dif- 
ference approximations.   A  total  of  nine  different  runs 
were  made,  each  of  which  applied  a  different  combination 
of  velocity  at  the  canard  exit  plane  and  ratio  of  mass  flow 
on  each  side  of  the  canard  ejector.   Each  combination  was 
obtained  by  altering  the  boundary  values  assigned  to  the 
canard  flaps. 

The  results  of  each  of  the  runs  were  similar  in  shape. 
Consequently,  only  one  of  the  outputs  is  included  subsequent 
to  the  program. 
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300 
400 
500 
900 


C 
C 

c 


10 


15 


20 


25 


30 


=  13,69,14 

N 

M 

1 

1 


FORMAT  (//) 

FORMAT  (6F8.2,6F4.2,//> 

FORMAT  (6F8.2f2F4.2,4Xf3F4.2»//J 

FORMAT  (/////} 

DIMENSION  A (4947) , B < 20 1 , 1 ) , C ( 20 1 } 

DIMENSION  DC201) 

INPUT  THE  ELEMENTS  OF  THE  COEFFICIENT  MATRIX 

DO  5  I  =  1,4869 

A( I )  =  0.0 

CONTINUE 

DO  10  I  =  1,201 

C( I)  =  0. 

CONTINUE 

M  =  2 

N  -  0 

DO  15  I 

M  =  M  + 

J  =  I  + 

K  =  J  - 

L  =  J  + 

A(K)  =  1. 

A(L)  =  1. 

A(J)  =  -10. 

K  =  J  +  12 

A(K)  =  4. 

N  =  N  4-  1 

CONTINUE 

M  =  2 

N  =  6 

DO  20  I 

M  =  M  + 

J  =  I  + 

K  =  J  - 

L  =  J  + 

A(K)  =  1. 

A(L)  =  1. 

A(J)  =  -4. 

K  =  J  +  12 

A(K)  =  1. 

N  =  N  +  1 

CONTINUE 

DO  25  I  =  223,348,25 

•  11 
12 
13 
24 

4. 

1. 

-10. 

1. 

4. 


=  112,168,14 

N 

M 

1 

1 


J  = 

K  = 

L  = 

M  = 

A(I  ) 

A(  JJ 

A(K) 

AIL) 

A  ( M ) 

CONTINUE 

DO  30  I  =  523, 

I  + 

I  + 

I  + 

I  + 


648,25 

11 

12 

13 

24 
4. 
1. 

-10. 
1. 
4. 
CONTINUE 
DO  35  I  =  823,948,25 

11 

12 

13 

24 
4. 
1. 


(I) 
(J) 
(K) 
(L) 
(M) 


J  =  I 
K  =  I 
L  =  I 
M  =  I 
A(I) 
A(J) 


i- 
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35 


40 


45 


50 


55 


60 


65 


DO  40 
J  =  I  + 
K  =  I  + 
L  =  I  + 
M  =  I  + 
A(I)  = 
A(J)  = 
A(K)  = 
A(L)  = 
A(MJ  = 
CONTINU 
DO  45  I 
J  =  I  + 
K  =  I  + 
L  =  I  + 
M  =  I  + 
A(  I)  = 
A(J)  = 
A(K)  = 
A<U  = 
A(M)  = 
CONTINU 


A(KI  =  -10. 
A(L)  =  1. 
A(M)  =  4. 

CONTINUE 

=  2248,2373.25 

11 

12 

13 

24 
4. 
1. 

-10. 
1. 
4. 
E 

=  3023,3123,25 

11 

12 

13 

24 
1. 
4. 

-10. 
4. 
1. 
E 

=  3323,3423,25 

11 

12 

13 

24 
1. 
4. 

-10. 
4. 
1. 
E 

=  3623,3723,25 

11 

12 

13 

24 
1. 
4. 

-10. 
4. 
1. 
E 

=  3923,4023,25 

11 

12 

13 

24 
1. 
4. 

-10. 
4. 
1. 
E 

=  4223,4323,25 

11 

12 

13 

24 
1. 
4. 

-10. 
4. 
1. 
E 

=  4523,^623,25 

11 


DO  50 
J  =  I 
K  =  I 
L  =  I 
M  =  I 
A  (  I  )  = 
A(J)  = 
A(K)  = 
A(L)  = 
A(M)  = 
CONTINU 
DO  55  I 


J  = 

K  = 
L  = 
M  = 
A(I) 
A(J) 
A(K) 
A(L)  = 
A(M)  = 
CONTINU 


DO  60 


J  = 

K  = 
L  = 
M  = 
A(I) 

A(J) 
A(K) 
A(L)  = 
A(M)  = 
CONTINU 
DO  65  I 


J  = 

K  = 
L  = 
M  = 
A(I) 
A(J)  = 
A(K)  = 

A«n  = 

A(M)  = 

CONTINU 
DO  70  I 
J  =  I  + 
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70 


75 


80 


85 


90 


95 


K  =  I  + 
L  =  I  + 
M  =  I  + 
A(I)  = 
A  (J)  = 
A(K)  = 
A(U  = 
A(M)  = 
C0N7INU 
DO  75  I 
J  =  I  + 
K  =  I  + 
L  =  I  + 
M  =  I  + 
A(I)  = 
A(J)  = 
A(K)  = 
A(L)  = 
ACM)  = 
CONTINU 


DO  80 


J  = 

K  = 
L  = 

M  ■ 
A(I) 
A(J)  = 
ACK)  = 
ACL)  = 
ACM)  = 
CONTINU 


DC 
J  = 

K  > 
L  ■■ 
M 
/8 
A 


85 
I 
I 
I 
I 


12 

13 

24 
1. 
4. 

-10. 
4. 
1. 
E 

=  398,498,25 

11 

12 

13 

24 
1. 
1. 
-4. 
1. 
1. 
E 

=  698,798,25 

11 

12 

13 

24 
1. 
1. 
-4. 
1. 
1. 
E 

=  2148,2223,25 

11 

12 

13 

24 
1. 
1. 
-4. 
1. 
1. 
E 


(I) 

(J) 
ACK)  = 
ACL)  = 

ACM)  = 

CONTINU 

DC  90  I  =  2848,2973,25 

+ 

-s- 


J  = 

K  = 

L  = 

M  = 

ACI) 

ACJ) 

ACK) 

ACL) 

ACM) 

CONT 


INU 


1 
1 
1 
2 

1. 

1. 

-4 

1. 

1. 

E 


DC  95 


J  = 

K  = 
L  = 
M  = 
ACI) 
ACJ)  = 
ACK)  = 
ACL)  = 
ACM)  = 
CONTINU 
DO  100 


=  3148,3273,25 


J  = 

K  = 
L  = 
M  = 
ACI) 
ACJ) 
ACK) 


1. 
1. 
-4. 
1. 

1. 

E 

I  =  3448,3573,25 
+  11 
*•  12 
+  13 
+  24 

:   1. 

'   1. 

■■    -4. 
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ACL)  =  1. 

ACM)  =  1. 
100  CONTINUE 

DO  105  I  =  3748,3873,25 

J  =  I  +  11 

K  =  I  +  12 

L  =  I  +  13 

M  =  I  +  24 

ACI)  =  1. 

A(J)  =  1. 

A(K)  =  -4. 

A(U  =  1. 

A(MJ  =  1. 
105  CONTINUE 

DO  110  J  =  4048,4173,25 

J  =  I  +  11 

K  =  I  +  12 

L  =  I  +  13 

M  =  I  +  24 

A(I)  =  1. 

A(J)  =  1. 

ACK)  =  -4. 

ACL)  =  1. 

ACM)  =  1. 
110  CONTINUE 

DO  115  I  =  4348,^473,25 

J  =  I  +  11 

K  -  I  +  12 

L  =  I  +  13 

M  =  I  +  24 

AC  I)  =  1. 

ACJ)  =  1. 

ACK)  =  -4. 

ACL)  =  1. 

ACM)  =  1. 
115  CONTINUE 

DO  120  I  =  2423,2523,25 

J  =  I  +  11 

K  =  I  +  12 

L  =  I  +  13 

M  =  I  +  24 

ACI)  =  1. 

ACJ)  =  1. 

ACK)  =  -4. 

ACL)  =  1. 

ACM)  =  1. 
120  CONTINUE 

DO  125  I  =  1048,1098,25 

J  =  I  +  11 

K  =  I  +  12 

L  =  I  +  13 

M  =  I  +  23 

ACI)  =  1. 

ACJ)  =  1. 

ACK)  =  -4. 

ACL)  =  1. 

ACM)  '  1. 
125  CONTINUE 

DO  130  I  =  1123, 1248,25 

J  =  I  +  11 

K  =  I  +  12 

L  =  1  +  13 

M  =  I  +  23 

ACI)  =  8. 

ACJ)  =  5. 

ACK)  =  -50. 

ACL)  =  5. 

ACM)  =  32. 
130  CONTINUE 

DO  135  I  =  1398, 1523,25 

J  =  I  +  1 

K  «  I-  +  11 
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L  =  I  + 
M  =  I  + 
N  =  I  + 
A(J)  = 
A(K)  = 
A(L)  = 
A(M)  = 
A(N)  = 

135  CONTINU 
DO  140 
J  =  I  + 
K  =  I  + 
L  =  I  + 
M  «  I  + 
N  =  I  + 
A(J)  = 
A(K)  = 
AIL)  = 
A(M)  = 
A(N)  = 

140  CONTINU 
A(1661) 
DO  145 
J  =  I  + 
K  =  I  + 
L  =  I  + 
M  =  I  + 
N  =  I  + 
A(J)  = 
A(K)  = 
A(L)  = 
A(M)  = 
A(N)  = 

145  CONTINU 
DO  150 
J  =  I  + 
K  =  I  + 
L  =  I  «- 
M  =  I  + 
N  =  I  + 
A(J)  = 
A(K)  = 
A(L)  = 
A(M)  = 
A(N)  = 

150  CONTINU 
DO  155 
J  =  I  + 
K  =  I  + 
L  =  I  + 
M  =  I  + 
A(I)  = 
A(J)  = 
A(K)  = 
A(L)  = 
A(M)  = 

155  CONTINU 
DO  160 
J  =  I  + 
K  =  I  + 

L  =  I  + 
M  =  I  + 
A(I)  = 
A(J)  = 
A(K)  = 
ACL)  = 
A(M)  = 
160  CONTINU 
N  =  5 
M  =  0 
DO  170 
L  =  I  + 


12 

13 

23 
16. 
1. 

-34. 
1. 
16. 
E 
I  =  1598,1648,25 

1 

11 

12 

13 

23 
4. 
1. 

-10. 
1. 
4. 
E 

=  0. 
I  =  1673,1798,25 


1 

2 

3 

24 


=  1948,2073,25 


0. 


=  2548,2673,25 


1 

2 

3 

24 


=  2723,2823,25 

1 

2 

3 


I  =  4648,4743, 19 
M 


80 


170 


175 


J  =  L  H 

•  11 

K  =  L  ^ 

■  12 

A(L)  = 

1. 

A(J)  = 

1. 

A(K)  = 

-4. 

J  =  L  -» 

13 

A(J)  = 

1. 

M  =  M  * 

•  N 

N  =  N  - 

•  1 

CGNTINUE 

N  =  3 

M  =  0 

DO  175 

I  =  47< 

L  =  I  4 

M 

J  =  L  + 

11 

K  =  L  + 

12 

A(L)  = 

1. 

A(J)  = 

4. 

A(K)  = 

-10. 

J  =  L  + 

13 

A(J)  = 

4. 

M  =  M  + 

N 

N  =  N  - 

1 

CONTINUE 

A(l)  = 

-10. 

A(2)  = 

1. 

A(13J  = 

4. 

A(99)  = 

1. 

A(100) 

=  -9. 

A(101) 

=  2. 

A(112) 

=  3. 

A(211J 

=  0. 

A (234) 

=  0. 

A(373) 

=  3. 

A(384) 

=  1.  , 

A(385) 

=  -9. 

A(386) 

=  2. 

A(397) 

=  3. 

A(511) 

=  0. 

A<534) 

=  0. 

A(673) 

=  3. 

A (684) 

=  1. 

A(685) 

=  -9. 

A (686) 

=  2. 

A(697) 

=  3. 

AC811J 

=  0. 

A(834) 

=  0. 

AC973) 

=  3. 

A(984) 

=  1. 

A(985) 

=  -9. 

A(986) 

=  2. 

A(997) 

=  3. 

A (998) 

=  1. 

A(1009) 

=  1. 

A(1010) 

=  -4. 

A(1011) 

=  1. 

AU022) 

=  1. 

A(1023) 

=  4. 

A(1034) 

=  3. 

A(1035) 

=  -18. 

A(1036) 

=  3. 

A(llll) 

=  0. 

A(1134) 

=  0. 

A(  1298) 

=  4. 

A(1309) 

=  3. 

AQ310) 

=  -18. 

A(1321 ) 

=  8. 

A(1273) 

=  8. 

AC1284J 

=  2. 

A(1285) 

=  -30. 

A( 1286) 

=  4. 

A(1296) 

=  16. 

4795,4837, 14 
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Al 

1324) 

- 

2. 

A! 

1335) 

=r 

-12. 

Al 

1336) 

= 

2. 

Al 

.1346) 

= 

4. 

Al 

1349) 

= 

4. 

Al 

il359) 

= 

3. 

Al 

1360) 

= 

-18. 

Al 

1361) 

= 

3. 

A< 

1371) 

= 

8. 

Al 

,1374) 

= 

4. 

Al 

1384) 

= 

3. 

Al 

1385) 

= 

-18. 

Al 

1396) 

= 

8. 

Al 

1409) 

= 

0. 

Al 

1549) 

= 

24. 

Al 

1559) 

= 

2. 

Al 

1560) 

= 

-54. 

Al 

1561) 

= 

4. 

Al 

1571) 

= 

24. 

Al 

1574) 

= 

3. 

A  I 

1584) 

= 

1. 

Al 

,1585) 

= 

-9. 

Al 

1596) 

= 

3. 

Al 

1609) 

= 

0. 

Al 

1661) 

= 

0. 

Al 

1684) 

= 

0. 

A< 

1824) 

= 

8. 

Al 

,1834) 

= 

4. 

Al 

1835) 

= 

-30. 

Al 

1836) 

= 

2. 

Al 

1846) 

= 

16. 

Al 

1849) 

= 

8. 

Al 

1859) 

= 

3. 

Al 

1860) 

= 

-18. 

Al 

1871) 

= 

4. 

Al 

1874) 

= 

4. 

Al 

1885) 

= 

-12. 

Al 

,1886) 

= 

2. 

Al 

,1897) 

= 

2. 

Al 

1899) 

= 

8. 

Al 

,1909) 

= 

3. 

Al 

1910) 

= 

-18. 

Al 

1911) 

= 

3. 

Al 

,1922) 

= 

4. 

Al 

1924) 

= 

8. 

Al 

[1934) 

= 

3. 

Al 

,1935) 

= 

-18. 

Al 

,1947) 

— 

4. 

Al 

11959) 

s 

0. 

Al 

2099) 

= 

3. 

Al 

'2109) 

= 

1. 

Al 

12110) 

= 

-9. 

A 

[2111) 

= 

2. 

Al 

,2122) 

= 

3. 

Al 

[2124) 

= 

1. 

Al 

[2134) 

r 

1. 

Al 

2135) 

= 

-4. 

A 

.2136) 

= 

1. 

A 

[2147) 

= 

1. 

A 

[2148) 

= 

0. 

A 

[2236) 

= 

0. 

A 

;2259) 

= 

0. 

A 

;2398) 

= 

3. 

A 

12409) 

s 

1. 

A 

[2410) 

= 

-9. 

A 

[2411) 

= 

2. 

A 

(2422) 

= 

3. 

A 

!2536) 

= 

0. 

A 

[2559) 

~ 

0. 

A 

(2698) 

- 

4. 

A 

(2709) 

- 

2. 

A 

(2710) 

- 

-12. 
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A(2711)  =  4. 

AC2722)  =  2. 

A(2836)  -  0. 

A(2859)  -  0. 

A(2998)  =  3. 

A(3009)  =  4. 

A(3010)  =  -18. 

A(3011)  =  8. 

A (3022)  =  3. 

A(3136)  =  0. 

A(3159)  =  0. 

A(3298)  =  3. 

A(3309)  =  4. 

£(3310)  =  -18. 

A(3311J  =  8. 

A(3322)  =  3. 

A(3436)  =  0. 

A(3459)  =  0. 

AC35S8)  =  3. 

A (3609)  =  4. 

A(3610)  =  -18. 

A(3611)  =  8. 

A(3622)  =  3. 

A(3736)  =  0. 

A(3759)  =  0. 

A(3898)  =  3. 

A(3909)  =  4. 

A(3910)  =  -18. 

A (3911)  =  8. 

A (3922)  =  3. 

A(4036)  =  0. 

A(4059)  =  0. 

A(4198)  =  3. 

A(4209)  =  4. 

A(4210)  =  -18. 

A (4211)  =  8. 

A(4222)  =  3. 

A(4336)  =  0. 

A(4359)  =  0. 

A(4498)  =  3. 

A(4509)  =  4. 

A (4510)  =  -18. 

A(4511)  =  8. 

A(4522)  =  3. 

A(4636)  =  0. 

A(4659)  =  0. 

A(4777)  =  3. 

A(4788)  =  4. 

A(4789)  =  -18. 

A(4790)  =  8. 

A(4857)  =  1. 

A(4868)  =  4. 

A(4869)  =  -10. 
C 

C      INPUT  THE  BOUNDARY  CONDITIONS 


C 


C(l)  =  -46. 
C(2)  =  -32. 
C(3)  =  -2  8. 
C(4)  =  -24. 
C(5)  =  -20. 
C(6)  =  -16. 
C(7J  =  -12. 
C(8)  =  -2.5 
C(9)  =  -2. 
C(10)  =  -1.5 
C(ll)  =  -1. 
C(12)  =  -.5 
C(13)  =  -10. 
C(25)  =  -10. 
C(37)  =  -10. 
C(45)  =  -16. 
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c 

c 
c 
c 
c 


c 
c 
c 
c 
c 


180 


210 


220 


230 


240 


C 

c 
c 
c 
c 
c 
c 

C( 
C( 

c 
c 

c 
c 
c 
c 
c 

C( 
C( 
C( 
C( 
C( 
C( 

c 

c 
c 
c 
c 
c 
c 

C( 
C( 

DO 
B( 

CO 


491 
56) 
57) 

60) 

67) 

68) 

78) 

79) 

90) 

71) 

82) 

94) 

106) 

118) 

130) 

142) 

154) 

166) 

178) 

190) 

191) 

192) 

193) 

194) 

195) 

196)  = 

197)  = 

198)  = 

199)  = 

200)  = 

201)  = 
180  I 

1,1)  = 
NTINUE 


-50. 

-6. 

-8. 

-10. 

-6.25 

-3.125 

-10.5 

-14. 

-3.5 

-50. 

-10. 

-10. 
=  -30. 
=  -10. 
=  -10. 
=  -10. 
=  -10. 
=  -10. 
=  -10. 
=  -19. 
=  -8. 
=  -7. 
=  -6. 
=  -5. 


-9. 
-2.5 
-2. 
-1*5 
-1. 
-0.5 
=  1,201 
C(I  ) 


SOLVE  THE  TWO  HUNDRED  AND  ONE  EQUATIONS 
SIMULTANEOUSLY  TO  OBTAIN  A  SOLUTION  FOR  THE 
STREAM  FUNCTION!  AT  DISTINCT  POINTS  IN  THE  FLOW 


001, IER) 


CALL  GELB(B, A, 201, 1,12,12, 
WRITE  (6,900) 
DO  210  I  =  1,201 
D(I)  =  B{ I  ,1) 

CONTINUE 


WRITE  THE  SOLUTION  TO  THE  STREAM 
EACH  DISTINCT  POINT  ALLOWING  FOR 
CONSIDERATIONS 


FUNCTION 
VARIABLE 


AT 
MESH 


I  =  1,37,12 


(D(L) ,L=J,K) 


11 


*  , 


L=J,K) 


DO  220 

J  =  I 

K  =  J  +  11 

WRITE(6,400) 

WRITE  (6,300) 

CONTINUE 

DO  230  I  =  49,71 , 

J  =  I 

K  =  J  +  10 

WRITE(6,500)  (D(L 

CONTINUE 

WRITE  (6,300) 

DO  240  I  =  82,94, 

J  =  I 

K  =  J  +  11 

WRITE(6,400)  (D(L),L=J,K) 

WRITE  (6,300) 

CONTINUE 

DO  250  I  =  106, 190,12 

J  =  I 

K  =  J  +  11 

WRITE  (6,400)  (D(L),L=J,K> 

WRITE  (6,900) 
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250  CONTINUE 
STOP 
END 


A  sample  output  follows  on  the  subsequent  page.   For 
this  case,  the  flow  velocity  was  evenly  distributed  in  the 
canard  augmentor  and  was  constrained  to  be  twice  the 
velocity  of  freestream  conditions. 

Boundary  conditions  are  not  included  in  the  output. 
The  decimal  point  of  each  stream  function  value  represents 
the  approximate  relative  position  of  each  solution  point  in 
the  domain  depicted  in  Figure  2. 
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